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The effect of fluctuations on the transport and thermodynamic properties of two-dimensional 
superconductors in a magnetic field is studied at low temperature T ^ T^o- The fluctuation conduc- 
tivity is calculated in the framework of the perturbation theory with the help of usual diagram tech- 
nique. It is shown that in the dirty case the Aslamazov-Larkin, Maki- Thomson and Density of States 
contributions are of the same order. At extremely low temperature T/Tco ^ {H — IIc2{0)) /IIc2{0) 
the total fluctuation correction to the normal conductivity is negative in the dirty limit and depends 
on the external magnetic field logarithmically Sa oc In {H — i/c2(0)). In the non-local clean limit, the 
Aslamazov-Larkin contribution to conductivity is evaluated with the aid of the Helfand-Werthamer 
theory. The longitudinal and Hall conductivities are found. The fiuctuating magnetization is calcu- 
lated in the one-loop and two-loop approximations. 

PACS numbers: 74.76.-w, 74.40. -f-k, 74.25.Fy 



I. INTRODUCTION 



Over the last decade there has been continuing inter- 
est in Quantum Phase Transitions. A particular atten- 
tion has been focused on two-dimensional systems which 
possess some unusual properties at low temperatures. It 
is remarkable, that a phase transition at zero tempera- 
ture is possible in the framework of the usual BCS theory 
of superconductivity. The transition temperature can be 
suppressed either by magnetic impurities or by magnetic 
field. It is interesting to find the fluctuation conductivity 
as a function of the closeness to the transition in these 
cases. The impurity-driven quantum phase transition has 
been considered by Ramazashvili and Coleman.EI Their 
consideration was based on the renormalization group 
analysis of the Aslamazov-Larkin correction to conduc- 
tivity. Fluctuations in an external magnetic field have 
been considered in different systems and various limit- 
ing cases.a u However, up till now, there is no consistent 
microscopic theory of superconducting fluctuations near 
Hc2{0)- The purpose of the present paper is to develop 
such a theory for two-dimensional superconductors in the 
dirty and clean limits. 

We begin with a brief review of studies of fluctua- 
tions in superconductors. The sulwect was initiated in 
the work of Aslamazov and Larkin.B The conductivity of 
fluctuating |Cooper pairs was calculated in zero magnetic 
field. MakiS and Thomsona included effects of electron 
scattering off the fiuctuations. It was found that there is 
another badly divergent contribution known as anoma- 
lous Maki-Thomson correction. Physically, this correc- 
tion is connected with the coherent scattering of the elec- 
trons by the impurities and analogous to the weak local- 
ization correction. The divergence can be removed by in- 
troducing a pair-breaking rate. Note, that experimental 
results at r ~ Tco can be described by the Aslamazov- 
Larkin term only. This sugestes that the pair-breaking 



rate is relatively large in real superconductors. Later, 
Thomson and Maki returned to the issue and evaluated 
fiuctuation correction io the normal conductivity in fi- 
nite fields. ThomsonQ evaluated paraconductivity for 
small fields T ^ Tco and large fields parallel to a twct 
dimensional superconducting sample. Ami and MakiEl 
considered a dirty three-dimensional superconductor put 
in an arbitrarily strong magnetic field having calculated 
the diagrams numerically. However, some technical sim- 
plifications that had been made in the paper (namely, the 
dynamic fluctuations had been neglected) make the re- 
sults inapplicable at very low temperature. Moreover, 
three-dimensional case is very different from the two- 
dimensional one, as shown in the present paper. Let 
us mention some r|elatively recent results in this field. In 
1993 Aronov et alB developed a theory of transport phe- 
nomena in the fluctuation region in the dirty, clean and 
superclean {lOcT 1) limits. Their consideration was 
based on the Ginzburg-Landau equations and, thus, is 
applicable for relatively small fields H ^ i?c2(0) only. 
Beloborodov et alH have calculated the fluctuating con- 
ductivity of a three-dimensional granular superconductor 
in the region close to Hc2{0). 

we con- 



Our paper is structured as follows. In Sec. II A 



sider a two-dimensional dirty sample Tcqt <^ 1 (where r 
is the scattering time). We calculate the total fluctua- 
tion correction to conductivity which is described by the 
standard set of diagrams (see Fig. 1). We derive an ana- 
lytical expression for the fluctuation conductivity in the 
region close to the transition line at low temperature i.e. 
at t = T/Teo < 1 and /i = (i/ - i?c2(T)) /ffc2(0) < 1. 
It is shown that in the case t ^ h the total correction is 
positive and has the usual form 5a oc Tcq [T — Tc{H)) ^, 
while at extremely low temperature t <^ h (at zero tem- 
perature, in particular) the total correction becomes neg- 
ative and logarithmically divergent Sa oc In h. 

In Sec. II B we address the issue of fluctuations in clean 
superconductors. This problem is more complex, since 
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the elements upon which the diagrams are built (current 
vertices, cooperons etc.) are non-local in the clean limit. 
We argue, that the corresponding operators can be found 
on the basis of the Helfand-Werthamer theory.El We apply 
this theory to our problem and calculate all the necessary 
values in the following limiting cases: Wc <C T or uJcT <C 1 
(where luc = eHc2{0)/m ~ Tco (Tco/sf) is the cyclotron 
frequency). This allows us to treat the magnetic field 
effects semiclassically. The curving of the classical tra- 
jectories is taken into account by comparison with the 
Drude conductivity. The longitudinal and Hall conduc- 
tivities are found. It is shown, that the fluctuation cor- 
rection to conductivity in the clean limit is similar to the 
one in the dirty limit, except for an additional cyclotron 
resonance-like pole of the second order which appears in 
the clean case. At the end of Sec. 11 B, we qualitatively 



discuss the effects of orbital quantization on the fluctu- 
ation conductivity, i.e. Shubnikov-de Haas oscillations 
which become essential at low temperatures T oj^. 

we calculate thermody- 



In Sections III A and IIIB 



namic properties of a superconductor. We find that mag- 
netization is logarithmically divergent in the first approx- 
imation and exceeds Landau diamagnetism. It is found, 
that in the clean case de Haas- van Alphen oscillations can 
become observable at high-enough temperature. Under 
certain circumstances the oscillating part of the fluctuat- 
ing magnetization represents the dominant effect. 

In Sec. IV we calculate the free energy and magnetiza- 
tion in the two-loop approximation for a dirty supercon- 
ductor. We find that the divergence becomes more severe 
in the higher orders in the perturbation theory. We dis- 
cuss the area of applicability of the results obtained. We 
find that the fluctuation region is determined by /i ^ Gi, 
where Gi ^ (si^r)^^, for low temperatures t <C /i but 
it becomes wider h < ^/G\t for relatively large tempera- 
tures t^h. 



II. FLUCTUATING CONDUCTIVITY 

A. Dirty superconductors 

The fluctuation correction to the conductivity beyond 
the Ginzburg region can be found in the perturbation 
theory. There are terms of three different types describ- 
ing the fluctuation conductivity in the first (one-loop) 
approximation. The first one is the Aslamazov-Larkin 
(AL) term (see Fig. 1.1) which is connected with the 
direct conductivity of the fluctuating Cooper pairs. The 
AL contribution to conductivity is positive. Since some 
fluctuating pairs appear above the transition, the number 
of normal electrons decreases. According to the Drude 
formula this leads to some decrease in the conductivity of 
the normal electrons. This contribution is known as the 
Density of States (DOS) term (see Figs. 1.5 and 1.8). It 
is clear, that this correction must be negative. The third 
term is the Maki- Thomson (MT) contribution (see Fig. 



1.2) which is connected with the coherent scattering of 
the normal electrons. The sign of the MT term is not 
prescribed. 

In the presence of impurities all these contributions 
must be averaged out over the impurities positions. This 
can be done in the Aamework of the diagram technique 
developed long ago.EiJ There is a standard set of diagrams 
to be considered in our problem (see Fig. pi). 








FIG. 1. Impurity averaging diagrams contributing to con- 
ductivity in the first (one-loop) approximation. 

These diagrams are built of the following elements: 
A solid line represents the one-electron Green function 
which in zero field has the form (in the momentum rep- 
resentation) : 
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(1) 



where £ — {2n+l)T:T is the fermion Matsubara frequency 
and = e(p) — is the one-particle excitation spec- 
trum. Here, we consider the quadratic spectrum. 

In the presence of magnetic field A(r) the Green func- 
tions change and contain the effects of orbital quantiza- 
tion. However, in the presence of strong disorder uJcT <C 1 
or at relatively high temperatures T ^ lOc, the discrete 
Landau levels are smeared out and the effects of magnetic 
field can be treated semiclassically. This means that the 
Green function in the coordinate representation can be 
written as 



Ge{vx, ra) =Si°)(ri -r2)exp l-ie^j A{s)ds\ , 



(2) 



where Q^^^ (r) is the Green function in zero field and the 
path of integration in Eq.(^ is a straight line. Let us 
note here, that the system of units h — c = ks = ^ 
is used throughout the paper. The magnetic field H is 
considered in the Landau gauge A = (0, —Hx). 
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Another element is the fluctuation propagator or inter- 
action in the Cooper channel (wavy line) . It is a diagonal 
operator in the Landau representation. The correspond- 
ing matrix element has the form:Ell 



NiO) 



47rr 



(3) 



where n corresponds to the n-th Landau level, (0) is the 
density of states per spin at the Fermi surface, Q — 2-KmT 
is the bosonic Matsubara frequency, corresponding to 
the total energy in the Cooper channel, VIh = AeDH 
and T> = ^v'^t is the diffusion coefficient. Note, that 
Eq.(^ is obtained from the expression for the fluctua- 
tion propagator in zero-fleld by the interchange of "Dq^ 
by Qh (n + 1/2), with q being the total momentum in 
the Cooper channel. 

The shaded vertices in the diagrams are Cooperons, 
which describe the coherent scattering of two particles 
off the impurities. i-The expression for this quantity has 
the following form:M 



Cn{£l, £2) 



(-£l£2) 



T ei - £2 



(4) 



where ei and £2 are fermion Matsubara frequencies, cor- 
responding to the electron energies. 

To calculate the total fluctuation correction to the 
DC conductivity we have to evaluate all the diagrams 
1-10 as functions of the external Matsubara frequency 
Lu = 2t:i>T, perform analytical continuation to the real 
frequency axis, take the limit lu and sum up all the 
contributions. The static term, corresponding to a; = 0, 
is cancelled out in the final result. 

In the vicinity of Tco (transition temperature in zero- 
field) only the AL and anomalous MT terms are impor- 
tant. The typical arguments are as follows. The point 
of superconducting transition is determined by the pole 
of the fiuctuation propagator (wavy line). The AL dia- 
gram contains two such lines. Thus, close to the transi- 
tion the corresponding contribution is the most singular 
one. Another singularity is due to the diffusion-like pole 
(~iuj + Vq^^ which appears in the MT term! (recall, 
that the MT process is connected with the coherent scat- 
tering of electrons). At small q and w ^ this yields a 
singular contribution. 

Another simplification which can be made at T ^ TcO 
is the possibility to neglect dynamic fiuctuations in the 
MT and DOS terms. This means that instead of evaluat- 
ing sum over the internal boson frequency D. we can just 
take the first term f2 = 0, which gives the most singular 
contribution. In the AL term the f2-dependence is con- 
sidered in the fiuctuation propagators only and neglected 
in the current vertices. 



The situation changes if a magnetic field is applied.i 
In this case, instead of integrating over q, we have to 
trace the corresponding operators over the Landau lev- 
els. The AL diagram contains only one singular fluctua- 
tion propagator £0 corresponding to the lowest Landau 
level, since the current vertex is not a diagonal opera- 
tor in the Landau representation. Moreover, the small 
terms Dq^ , which exist in zero-fleld, have to be replaced 
by {n + 1/2) ^ Tcq. Obviously, the anomalous MT 
term does not possess any additional singularity in this 
case. Thus, we conclude, that different diagrams should 
give contributions of the same order if a large magnetic 
fleld is applied. 

Let us now perform a representative calculation on the 
example of the AL term (see Fig 1.1). The corresponding 
expression for the longitudinal component of the electro- 
magnetic response tensor has the following form: 



n=0 



£„(r2)£„+i(ri - uj) + £„($! - a;)£„+i(ri) 



(5) 



where factor 4 is due to the spin, constant c = 
AttN{0)'Dt'^ appears as a result of the integration over 
^ in the local current vertex (see (|A9|)), — eH/n is the 
number of states per unit area of a full Landau level and 
7r„„+i = {n\{-i\7 + 2eA{r)) J n + 1) = y^{n+l)eH are 
matrix elements of the kinetic momentum, u is the Mat- 
subara frequency corresponding to the frequency of the 
external electric field, and e are the internal bosonic 
and fermionic Matsubara frequencies respectively. 

As we have already mentioned, the main singularity 
comes from the fluctuation propagator corresponding to 
the lowest Landau level. Close to the transition it can be 
written as 



Com 



1 



1 



A^(o) h + 2\n\/nH' 



(6) 



where h = (H - Hc2{T)) /Hc2{0). Let us note, that 



2tt 



TcO a-nd the bosonic frequency 



n is of the order of temperature. Thus, we conclude that 
at very low temperatures t <^ h we can replace the sum 
over in Eq.(^) by an integral. At relatively high tem- 
peratures t ^ h we can keep the first term in the sum 
only. If i ~ /i we have to evaluate the sum. This also 
means that we have to consider the effects of dynamic 
fluctuations as well. 

Let us discuss some simpliflcations that can be made 
in our case {t <C 1). First of all, we can consider only 
the first term n = in the sum over Landau levels in 
Eq. (^ . Only this term give a singular contribution com- 
ing from Co- Next, we see that the sum over the internal 
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frequency in (||) is determined hy T fin- This al- 
lows us to make expansions with respect to fl/riH ~ t ev- 
erywhere except £o with being a Matsubara frequency, 
. With the same accuracy, we can replace the sum over 
the fcrmion energy e in Eq.(||) by an integral. 

Evaluating the integral over e, we obtain from Eqs.(^, 

I) 



3nH/2 + \n-uj\ 



1 



In 



VLh-uj \VLHl2+\n-Lj\ + \uj\ 



1 



f2jy -|- UJ 



In 



3^7g/2 + 1»| + \uj\ 
nHl2 + \n\ 



(7) 



Now, we have to perform analytical continuation in the 
expression for the current response operator (||). In do- 
ing this, we can present the sum over the Matsubara 
frequency as an integral over the real frequency with the 
function coth which is chosen to generate poles at 

the points 2nimT\}3 Making use of Eqs.(^ - ^ we obtain 
the following expression for the conductivity (within the 
logarithmic accuracy) 



Sai = lim 



Qi{u! iuj) 



iui^O —IUJ 



ai I dn coth — 



n 



^2 + (11^)2 -^^^J 2T sinh^ # + (il^) 



dn 1 



, (8) 



ai = 


4 
3' 




Pi = 


2; 


a2 = 


-2, 




/32 = 


2; 


O-z = 


an — 


2 

3' 


/33 = 


/34 = 0; 


"5 = 


"8 = 


3 

2 ' 


/?5 = 




o-e = 


ag = 


5 

3 ' 


/36 = 




ai = 


"10 = 


1 

2 ' 


/37 = 




atot = 


2 

3' 




Aot = 


_ 8 
3' 



(12) 



where indexes correspond to a diagram number in Fig. 
1 and atotj Aot describe the total correction to the con- 
ductivity, which can be written as: 



5a = 



2f_ 



r 3 

-\n - - — + ip (r) + A{rtp' {r) - 1) 
n Zr 



(13) 



Let us consider some limiting cases. If the temperature 
is relatively large t ^ h, we obtain the following formula 
for the fluctuation conductivity: 



276^ t 



(14) 



If i7 < Hc2{0), we can introduce Tc{H) and rewrite 
Eq.(|l4[) in the usual way 



where ai = 4/3 and /3i — 2 are just numbers. One 
can see that the first integral in Eq. (j|) is logarithmically 
divergent. This divergence appears as a result of our ex- 
pansions in t. Thus, it has to be cut off at fimax ~ Tcq. 
The integrals in Eq.(H) can be easily calculated. The 
result is 



3e2 



T, 



cO 



2jTT'^h T-TciH)' 



(15) 



li H > Hc2{0), in the low-temperature limit t ^ h we 
have 



5a 



ala{h,t) + pifsiKt) 



(9) 



Sa 



2e^ 1 



(16) 



with 



and 



/„(/i,<)=lnJ--^-7^(r) 
n 2r 



l0{h,t)^r^' (r) 



1 

2^ 



(10) 



(11) 



where ^ — j and 7 — 0.577 is the Euler's constant. 

The other diagrams can be calculated analogously. 
The corresponding contributions to conductivity can be 
written in the same form as Eqs.(|8| — 11). Below we give 
the results in terms of constants a and f3: 



One can see, that even at zero temperature a logarithmic 
singularity remains and the corresponding correction is 
negative. 

Let us note, that the fluctuating conductivity depends 
on the magnetic field and temperature via their ratio 
h/t. The behavior of the conductivity in the vicinity of 
the critical point H = i?c2(0), T = depends on the 
way how one approaches this point. If the transition 
is driven by the magnetic field and the temperature is 
zero, than the fluctuating correction is negative and log- 
arithmically divergent. If the magnetic field is fixed and 
H < Hc2{0), than the correction is positive and diverges 
as {T — Tc{H)) ^ . In the other cases, there is a crossover 
between these two regimes. 
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0.025 0.05 0.075 0.1 0.1250.150.175 0.2 

Magnetic field h{0) = {H - H^iiO)) /H^iO) 
FIG. 2. Fluctuating conductivity ( p^ as a function of 
magnetic field is plotted for four different temperatures. 
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0.2 0.3 0.4 0.5 

Temperature t = T/Tco 
FIG. 3. Fluctuating conductivity ( p^ as a function of 
temperature is plotted for five different magnetic fields. 

The magnetic field dependence of the fluctuating con- 
ductivity is presented on Fig. 2. One can see, that if the 
magnetic field is relatively large, than the total correc- 
tion is negative. For any finite temperature, there is a 
region close to i?c2(0) where the correction is positive. 

The temperature dependence of the fluctuating con- 
ductivity is shown on Fig. 3. It is interesting, that the 
conductivity is a non-monotonous function of tempera- 
ture if the magnetic field exceeds Hc2{0)- 



B. Clean limit 

In this section, we investigate the fluctuation correc- 
tion to the conductivity in the limit TcT 3> 1. In this case, 
usual expressions for the particle-particle bubble, fluctu- 
ation propagator and current vertices are inapplicable. 
To calculate the diagrams we have to find these quan- 
tities in the presence of the magnetic field while taking 
into account their non-local structure. There are several 
effects associated with the magnetic field applied. First 
of all, the superconducting transition itself is governed 
by the magnetic field at low temperature. Another ef- 
fect is Shubnikov-de Haas oscillations in the conductivity 
due to the quantization of the energy levels. However, if 



LUcT <C 1 or T ^ cjc the oscillating terms are exponen- 
tially small and can be neglected. Note, that 



ei/c2(0) 



Tco 



Tco 



(17) 



In our formal derivation, we assume that either <^ 1 
or LUc T <^ TcQ. This allows us to consider low temper- 
atures without dealing with de Haas oscillations in the 
Green functions. The effect of the orbital quantization 
on the fluctuation conductivity will be briefly discussed 
at the end of this Section. Moreover, there is a purely 
classical effect due to the Lorentz force acting on the 
electrons forming fluctuating pairs. Namely, the mag- 
netic field results in curving of the classical trajectories. 
This curving leads to the cyclotron resonance and Hall 
effect in the fluctuation conductivity. First, we consider 
fluctuations neglecting the curving, which is eligible if 
ujcT ^ 1- Using the result obtained, we will be able to 
derive the formula valid in the superclean case locT ^ 1 
as well. 

We now proceed to calculate different blocks in the 
diagrams. Our calculation is based on the well-known 
Helfand-Werthamer theory developed long ago. In the 
seminal paper,Ll Helfand and Werthamer evaluated the 
matrix element Co for the Cooperon in a magnetic field, 
which determines the upper critical field Hc2{T). They 
proved the following mathematical statement which we 
will refer to as the Helfand-Werthamer (HW) theorem 
throughout the paper. 

Let us consider an operator O. Suppose, that its kernel 
in the coordinate representation has the following form: 



0(r, r') = 0(r - r') exp \^2ie J A(s)ds | . (18) 
Than, the operator can be written as 



0= / 0{r)e- 



(19) 



where tt — {p — 2ieA(f)) is the kinetic momentum, 
which can be expressed in terms of the creation and an- 
nihilation operators in the Landau representation, and D 
is the dimensionality of the system (13 = 2 in our case) . 

One can see that all the operators involved in our cal- 
culations satisfy the HW theorem. Namely, the particle- 
particle bubble current vertex 70(17,^) and the 
four Green function blocks i?Q^(r2,cj) in the coordinate 
representation can be written as a product of a function 
of the coordinate difference and the gauge factor. In the 
temperature range under consideration, we can treat the 
magnetic field effects semiclassically which means that 
the first factor O in Eq. (|l8|) can be considered in zero 
field. 

To calculate the matrix elements of our interest we 
will do the following. First, we calculate an operator in 
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zero-field in the momentum representation 0(q). We ap- 
ply the Fourier transformation to this function and put 
the value obtained 0(r) in Eg. (|l9|) . Than, we evaluate 
the matrix elements for this operator and perform the 
integration over r. Finally, we perform the frequency 
summation left (over the fermion energy e). 

Let us start with the calculation of the non-local fluc- 
tuation propagator which has the form: 



c{n) = 



1 



where g is the interaction constant and 



fl(q,f})-T^n,(q,r!), 



(20) 



(21) 



wit h the particle-particle bubble He (q, il) defined by 
Eq.(All). Note, that in the clean limit we can neglect 
the impurity dependence in n£(r2) and in the fluctuation 
propagator. 

The matrix elements can be calculated by expressing 
TT in terms of the creation and annihilation operators 
and a and expanding the exponentials.113 One obtains 



exp (— ir-Tr) 



eL V- {-ip) 

^ k\l\ 

k. 1=0 



(22) 



where p — r/^/2rH and — y/2eH is the magnetic 
length. Due to the integration over (f> only the diago- 
nal matrix elements survive and we have the following 
expression 

oo 

^ 

(23) 

where L„ is the Laguerre polynomial of the n-th order. 

At low temperature, we can replace the sum over s by 
an integral and we have 



n„(r!) = iv(o) 



In (2A) 



(-1)" f dx In (a + + e^"^L„(2x) 



(24) 



where we have introduced A = ilit^ beine; the lower 
Hmit of integration over e and A = 2niitiD_ -^j^^ich is the 
BCS high-energy cut-off. Obviously, A ~ T/Tco <C 1 and 
A> 1. 

Let us realize, that to find the most singular contri- 
bution to the conductivity we need to know no(ri) and 



Hi (Jl) only. Making expansions with respect to A in 
Eq. (p4[) , one obtains 



and 



no(f^) = iV(0) [In (2/7A) - V^A] (25) 



Hi {9) = N{{)) [1 + In (2 V7A) + A^] . (26) 



Thus, the fluctuation propagator, corresponding to the 
lowest Landau level can be written in the vicinity of the 
transition as follows 



Com = 



1 



N{Q)h+^\n\/T,^' 



(27) 



where T^o is the transition temperature in zero field, 

h = {H-H,2{T))/H,2{0) and eH.^iO) = ^ 
is the upper critical field at zero temperature. 

The current vertex can be evaluated in the same fash- 
ion as the fluctuation propagator. However, the cor- 
responding calculation is more cumbersome (see Ap- 
pendix) . In the momentum representation the vertex has 
the form: 

7(q; n, oj, e) = 4 ^■'(^' ^)/(^- <?) ' (28) 

^ i=l 

The corresponding functions an d co n stan ts are defined 
in the Appendix by formulae (|A2|— |A8|). Eq.dl) is 
large mostly due to the thcta-functions. However its q- 
dependent part has a simple form ^f{q), where /(O) = 
and /(oo) is finite. 

Making use of the Helfand-Werthamer theorem, we ob- 
tain the following expression for the operator 7: 



i=l 



dqf{et,q) 



27T 



f f d(hY 

X / drrJi(qr) / exp (— ir-fr) (29) 

J J 2TTr 



Evaluating matrix elements and performing integration 
over and r, one can see, that the current vertex pos- 
sesses only near-diagonal non-zero matrix elements and 
they have the form: 



{n\r{n,u,e)\n+l) 



n + 1 



4 

E^^^^^^/'^'^V(e..)e~^«'^L«(24,^). (30) 

1=1 n 



In order to calculate the most singular contribution to 
the conductivity we have to know 701 only. Taking the 
corresponding integral by parts one obtains 



6 



701 



(0,w,e) = V2TTNiO)rHj2 



and (Q) together and using the following formula for the 
current response operator 



dq 



(31) 



To evaluate the remaining integrals we have to use the ex- 
plicit expressions for the function s /( e, q) and rii{e,Q,uj) 
(see Appendix). Using formula (|A2|), one obtains after 
elementary integration over q 



70 



1 (rj, uj, £ 



i=l 



1- V^|(5i|e*'erfc|(5, 



(32) 



where we have introduced 5i — tirnjvp. Now, we have 
to perform the summation over the fermion frequency e. 
For T <C Tco we can replace this sum by an integral over 
(5 taken in the a ppro p riat e limits well-defined by the 0- 
functions in Eqs.(|A^ — A8). The corresponding indefinite 



integral can be easily evaluated and we finally obtain 



7oi(ri, w) 



1 

2^ 



de7oi(r2,Ci^,£) 



> — sgn di e • erfc \di 

1=0 



limits ; 



(33) 



where the limits of integration are determined by the 
theta- functions in r]i. Let us note, that at 5 = c» the 
expression written above vanishes, while the other limits 
of integration are such that S ~ T/Tco <C 1. 



There are four terms in Eq.(A3). The last two contain 
the factor cj~^sgn£ , while the first two are proportional 
to the following factor 



1 + ujT sgne 



sgne 



1 - 



1 



1 + LOT Sgne 



Using Eq.(p3[), one can see, that the tj^^-terms are can- 
celled out exactly. Thus, the current vertex can be writ- 
ten (we keep only linear terms with respect to the fre- 
quencies) : 



7oi(f^,t^) = - 



iV(0)rg 
V2 



1 



\/7r rn 

2 VF 



1-1- |w|r 

{\n\ + \n-uj\ + \uj\) + o{t) 



(34) 



We see, that the current vertex is proportional to 
(1 -|- Ic^It) ^ and its frequency dependence is determined 
by the two pairs of 0-functions in Eq.(Al). These terms 
exist, when the poles e and (e— w) are located in the oppo- 
site half-planes of the complex plane Let us note, that 
the similar situation takes place when calculating Prude 
conductivity of the normal metal. Getting Eqs.(E^), (|2£ 



Q{u;) - Sz/e^T V7o\(r!,c^)/:o(f^)A, 



(35) 



we can calculate the AL contribution to the conductivity. 
Note, that in the framework of our approximation, we can 
treat £1 as a constant, since it does not have any linear 
ri-dependence (see Eq.(|2^)). The analytical continuation 
yields the following expression for the conductivity (valid 
within the logarithmic accuracy): 



Sa 



5(7 



1 



(1 — iujr)'^ (1 — iujr)'^ 



I^(h,t)+Ip{h,t) 



(36) 



where functions Ia{h,t) and If}{h,t) are defined by 

Eqs.(p^, |ll| ) with parameter r ~ \J^l^ which differs 

by a constant from the one in the dirty case. 

Equation ( p6| ) is valid for uj <C T^o. Recall, that the 
Drude conductivity has the following form: 



Co = 



1 



1 — iujT m 1 — iujT 
Thus, the total longitudinal resistivity reads 



6a 



(37) 



(38) 



We see, that the fluctuation correction to the resistiv- 
ity does not depend on the external frequency (unless 
Tco) and can be considered as a correction to the 
collision integral i. Physically, this means, that the AC 
electric field acts on the normal electrons, rather than 
superconducting fiuctuations. 

In the superclean case uJcT < 1, we have to take into 
account the curving of the classical trajectories. This 
curving results in the Hall term in the conductivity and 
cyclotron resonance-like effects. The Hall term can be 
written as 

_ m 

Pxy — :t + "Pxyj 

where the second term is due to the superconducting fluc- 
tuations. The reasonings described above suggest that 
this term, which describes the curving of the fluctuat- 
ing pairs, is of the order of lOc/Tco and can be neglected. 
Hence, calculating inverse matrix we find the follow- 
ing formula for the fiuctuation conductivity 



1 



I — i (uj ± LUc) t] 



rSa, (39) 



where Sa is defined by Eq.(p6|). Let us note, that Sa rep- 
resents the longitudinal conductivity with no respect to 
the curving. The corresponding Hall term Saxy can only 
appear in the presence of a particle-hole asymmetry. It 
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does not exist in the framework of our approximation. In 
the paper of Aronov et aln this additional Hall term was 
controlled by the phenomenological parameter Tc ^g" J'' ■ 
Let us now discuss the contributions coming from the 
MT and DOS diagrams. The electromagnetic response 
tensor can be written in the following form: 



(40) 



SI, 6 



where Ba^ represents a four Green function block. Let 
us consider this quantity on the example of the MT term. 
In the coordinate representation within semiclassical ap- 
proximation it has the form: 



Ba/3 (e, ^; r, r') = Bap (e, w, fi; r - r') 

X exp I —2ie J A(s)ds 



(41) 



where 

Bap (e,w,fi;r) 



Varp f dpdk 



J (27r)^ 





Ji{pr)Ji{kr) 



X Ge{p)Ge-u{p)Gn-e{k)gn+^-e{k). (42) 

Putting this expression into the HW theorem ( [l9| ) , evalu- 
ating the diagonal matrix element for n = and perform- 
ing integration over r we obtain the following expression 



(0|B,, {e,uj,n)\Q) = - 



dpdk 2,2 
27r ^ 



X Qe ip)Ge-uj {p)Qn-e (fc)5o+u;-e (fc)^O {P, k) , 

where 

^r%pk 



looip, k) = exp [-rl {p^ + e)]h 



(43) 



(44) 



and /i is the modified Bessel function of the first or- 
der. Since p ^ k ^ pp and, thus, r|^pfc ~ s/oJc ^ 
1, we can take the asymtpotical form of the Bessel 
function and then approximate the resulting exponent 



exp 



-rlip-ky 



by the (5-function. Thus, we have 



/oo(p, k) « -5[p- k). 
P 



(45) 



Performing integration with respect to k and introducing 
the density of states at the Fermi-surface A^(0) we obtain 
the usual expression for the 4-Green function block 



( \B.,^{e,u,9)\Q) = ^\lN{Q) 



j d£^pQe{p)Qe-u{p)Qn-e{p)Q'. 



(46) 



One can see, that the derivation of this expression does 
not depend on the purity of a superconductor. It is valid 
for the dirty and clean limits and any magnetic fields 
applied, unless uJc ~ sp- 

Let us note that Eq.(|4^) for the MT diagram and the 
similar expressions for the DOS diagrams are identical to 
the ones in the vicinity of T^q and do not involve mag- 
netic field at all. It is known, that DOS andJMT terms 
are strongly compensated in the clean limittll and this 
compensation takes place at the level of Green functions 
[i.e. before integrals over q). This suggests that, in the 
clean limit the only remaining diagram is the AL term 
even in the case of strong magnetic field. 

Let us now discuss quantum oscillations in the fiuctua- 
tion conductivity. At very low temperature these oscilla- 
tions become important. In this case, each current vertex 
contains an oscillating part. This oscillating part can be 
found by comparison with the Drude conductivity and 
can be written as 7 = 70 +7osc, where 7osc/7o ctosc/cto 
with CTosc being the oscillating part of the normal conduc- 
tivity (see e.g. Ando et altB) However, there are other 
"sources" of quantum oscillatiops. .The transition tem- 
perature Tc{B) oscillates as welOEZi and affects the fluc- 
tuation conductivity. Let us also realize, that the oscilla- 
tions of magnetization (de Haas- van Alphen oscillations) 
can influence Shubnikov-de Haas oscillations. Under cer- 
tain conditions, this effect may be dominant. Moreover, 
magnetization fluctuates as well and in the vicinity of 
the transition the fluctuations can exceed Landau dia- 
magnetism (see section 3 and Ref. ^8|). We see, that 
the oscillating part of the fluctuating conductivity has 
a complicated structure and can differ significantly from 
the usual Shubnikov-de Haas oscillations. 



III. THERMODYNAMICS. FLUCTUATING 
MAGNETIZATION 



A. Dirty case 

Considering thermodynamic properties of a film, we 
can calculate the free energy directly. In the cpe-loop 
approximation, the free energy can be written agl3 



Fi--r^Trln(l-5C'(17)), 
o 



(47) 



where (7(17) is the cooperon. 

Using Eqs.(||), (^ and (47), one can easily obtain the 
magnetization 

where d is the thickness of the film or the interlayer dis- 
tance, V = eB I'K is the number of states of a Landau 
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level and function Ia{h,t) is defined in Eq.(|lO[). Thus, at 
low temperature t <ti h the susceptibility takes the form: 



Xi = 



dH TT^hc'^ d 

One can see, that the fluctuation susceptibility (^) is 
large compared to the magnetic susceptibility of the nor- 
mal metal x l even far from the transition: 



Xi 



Gih 



(50) 



where Gi = (eft) ^ is the Ginzburg parameter. 



B. Clean case 

The calculation of magnetization in the clean limit can 
be done in the same fashion as in the dirty limit. How- 
ever, there are some features specific for the clean case. 
As we have already mentioned, de Haas oscillations be- 
come essential at low temperature in pure samples. These 
quantum oscillations appear in all quantities including 
Green functions, transition temperature Tc{H), fluctu- 
ating conductivity etc. The oscillating terms are pro- 
portional to the factors exp ^— and exp ^— 

Hence, the oscillations are strongly suppressed, unless 
ujcT ^ 1 and T/u!c ^ 1. Let us note, that de Haas- 
van Alphen oscillations in the magnetization can reveal 
themselves much earlier than the quantum oscillations in 
the other quantities. This is because magnetization is a 
derivative of the free energy with respect to the magnetic 
field. Even though the oscillating terms in the Green 
functions are small, they contain fast-oscillating func- 
tions cos {2tt£ p/ijjc) which may lead to observable effects 
in the oscillating magnetization. It is easy to calculate 
the fluctuating magnetization with respect to these ef- 
fects. 

We can use Eq. (E^) 



Fi = Tr In (l - gn{VL) 



(51) 



where H = Hq + Hosc is the particle-particle bubble which 
contains an oscillating part. The matrix element for the 
monotonous part of the particle-particle bubble corre- 
sponding to the lowest Landau level was found in Sec. 
II B (see Eq.(25)). The oscillating part has beejxxonsid- 
ered in a number of papers and has the formil3'E3 



Ho 



-87r3/27V(0) 



T 



cos I 2tt — 



, pLeH\ ( A 
X cos I 67r exp 



(52) 



where A = Gtt (ttT + ^) and /ie is the magnetic moment 
of an electron. For the sake of simplicity, we keep the first 
oscillating term only. 



In the vicinity of the transition we can present the 
magnetization in the following way 



(49) Afi = 



1 



d 



Tv 



H,2io) ^ Coin)-' -'nose dh 



[Coin)-' - Ho 



(53) 



where for Co(^Y) see Eq.(p7|). 

From Eqs.(p7|), ( ^7|) and (|5^), we obtain the following 
expression for the fluctuating magnetization 



Ml 



1 TrOl' 



In th 

t ^ 



1 h 



1 + 32^5/2 Zj^3in (^^e_F_ 



UJr 



X COS I DTT I exp 



A 

LOr. 



(54) 



3/2 



Let us note, that if T ^ uJc^ than T ^fsp jujc 
ep/Tci) S> 1 and the numerical factor in the oscillating 
term in Eq.([5^) is very large. Thus, we conclude that 
de Haas-van Alphen oscillations in magnetization may 
exist even in the absence of the Shubnikov-de Haas oscil- 
lations and oscillations of the transition temperature. It 
is worth mentioning, that the fluctuation effects exceed 
Landau diamagnetism in the clean limit as well (formula 
( pO| ) is valid with Gi = ef/Tcq). Thus, under certain 
circumstances (A ~ 1) the oscillating part of the fluctu- 
ating magnetization may be more important than the 
monotonous part of magnetization and the oscillating 
part in the Landau term. 



IV. TWO-LOOP APPROXIMATION. 
APPLICABILITY OF THE RESULTS. 

In the previous sections we found the fluctuation cor- 
rection to the the transport and thermodynamic proper- 
ties of a superconductor in a magnetic field in the first 
(one-loop) approximation. The purpose of the given sec- 
tion is to find the order of the subleading corrections. 
This will determine the area of applicability of the re- 
sults obtained. We shall calculate the magnetization in 
the two-loop approximation for a dirty superconductor. 
This correction can be easily calculated in view of the 
simplifications described above. 

In the two-loop approximation we have to deal with 
diagrams presented on Fig. 4. The corresponding con- 
tribution can be written in the coordinate representation 
in the following way 



xis:£(ri, r2; rs, r4)£n(ri, r2)>Co' (rs, r4), (55) 
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where Kg is the operator corresponding to the square 
blocks in the diagrams presented on Fig. 4. This oper- 
ator is famiUar froci the usual BCS Jtheory. It has been 
calculated by Makic^ and Caroli et aLtil and has the form: 



if£(ri,r2;r3,r4) 



7r7V(0) 



X 5{yi - Y2)5{yi ~ r3)(5(ri ^ r4) 



..1 .-I + I^^W 



-9(3)]% [9(2) 
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where we make use of the Maki's notations: 



(56) 




FIG. 4. Diagrams contributing to the free energy in 
the two-loop approximation. Similar diagrams appear in the 
derivation of the Ginzburg-Landau equations from the micro- 
scopic theory. 



5, 



\k) 



-iV- 2e(-l)'=A(r). 



In the coordinate representation, the fluctuation propa- 
gator can be expanded on the basis of the eigenfunctions 
in the magnetic field and has the form: 



/:n(r,r')= j ^X]'C"(")^np«W^"f«(^')' (57) 



where £„(f2) are matrix elements of the fluctuation prop- 
agator in the magnetic field (see Eq.(||)), ■!/)„p^(r') is the 



eigenfunction for an electron in a magnetic field in the 
Landau gauge and Py is the y-component of the momen- 
tum, which determines the orbit's center. Again, in the 
vicinity of the transition line we keep the n = term 
only in Eq. (|57|) . From Eqs.(|5| — we obtain the free 
energy per unit volume 



V 



7rA(0) 
2d 



1 



{\e\^^lH/2y 



Thus, the magnetization takes the form: 

M 1 dlliKt) 

^ (27r)2dA(0) i7c2(0) dh ' 

At low temperatures t <^ h we have 

^2 1 1 1 



(58) 



(59) 



(60) 



We see, that the second order correction is negative. 
From Eqs.(E^) and (|59|) we obtain the ratio 



M2 

Ml 



Gi 
27r 



t 



1 



2-ft 



1 h 

2^1 



(61) 



where Gi is the Ginzburg parameter. The one-loop ap- 
proximation is valid unless this ratio becomes of the order 
of unity. At low temperatures t <^ h, Eq.(|6l|) yields the 
following condition 



U t ^ h, we have 



/i> Gi. 



h > 



(62) 



(63) 



This indicates that at large enough temperatures the fluc- 
tuation region becomes wider. 

These results stand for the kinetic coefficients as well. 
In the clean case the formulae (|6^) and ( |63| ) are valid 
with Gi ~ ep/Tco- However, explicit calculations are 
more complicated due to the non-local structure of the 
K-operator. 

Let us note, that at an exponentially low temperature 
some other effects may reveal themselves. In the_dirty 
case, mesoscopic fluctuations may be important.E3 Re- 
ally, the upper critical field depends on disorder. The 
distribution of impurities is random. There are some re- 
gions where the concentration of the impurities is such 
that the upper critical field is smaller than the bulk value. 
These regions may form superconducting islands weakly 
coupled one with another. At extremely low tempera- 
ture the proximity effect and the Josephson coupling can 
make these mesoscopic fluctuations observable. The ef- 
fects due to the mesoscopic fluctuations will be consid- 
ered elsewhere. 
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V. CONCLUSION 



APPENDIX A: CURRENT VERTEX 



The central result of the paper is the existence of the 
logarithmic correction to the conductivity which persists 
down to zero temperature. This correction is shown to be 
negative in the dirty case. The minus sign comes from 
the DOS diagrams as well as from the anomalous MT 
term. The AL contribution is positive but numerically 
smaller. Let us note, that similar results (negative fluc- 
tuation correction to conductb/ity) exist for the granular 
and layered superconductors.OO In these cases the AL 
and MT contributions are parametrically small compared 
to the DOS term. 

The fluctuating magnetization exceeds conventional 
Landau diamagnetism for a very large range of fields. It 
is shown to be logarithmically divergent as well at T — > 0. 

Let us note, that singular behavior of transport and 
thermodynamic quantities at low temperature is due to 
the low dimensionality of the system. In the three di- 
mensional case the leading correction to conductivity is 
not singular Jctsd ^ V^- 

The results obtained in the present paper can be 
checked experimentally by measuring fluctuation conduc- 
tivity in two-dimensional and quasi-two-dimensional sys- 
tems. The results obtained in the dirty limit can be 
checked by measuring magnetoresistance in dirty super- 
conducting films at low temperature. In this case, there 
could be some experimental difficulties connected with 
the Hc3-eSects that can screen the bulk properties of a 
film. The edge effects can be excluded, for example, by 
putting a sufficient amount of magnetic impurities on the 
edge of the film. 

The clean case may be relevant to high-Tc supercon- 
ductoraS^ and, probably, to the recently discovered two- 
dimensional organic superconductors .t3 Let us note, that 
our results assume s-pairing and isotropic Fermi-surface 
which is not true for High-Tc superconductors. However, 
it can be shown, that the logarithmic singularity In ^ re- 
mains for any pairing type (with a coefficient different 
from our case). It is worth mentioning, that in the over- 
doped High-Tc superconductors the Ginzburg parameter 
Gi is small and, thus, the fluctuations are negligible. In 
the underdoped superconductors the fluctuations are ex- 
tremely large and they lead to a large pseudogap which 
makes the conventional Fermi-liquid theory inapplicable. 
Hence, optimally doped superconductors should be used 
to check the results obtained. 
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In this Appendix we derive the formula for the cur- 
rent vertex depending on two frequencies O and u; in the 
non-local clean limit. The corresponding result is used 
when calculating the Asla maz ov-Larkin contribution to 
the conductivity (see Sec. |II B| ). 

The current vertex is the triangle block in the AL dia- 
gram (see Fig. 1.1). It consists of three Green functions. 
In the momentum representation, it can be written as 

f (Pp 

7e(q;f7,cj) = J ^^-^vG£(p)Ge_^(p)Gn_e(q-p) 
= -iV(0) / d^- 



1 



V 



I [e ~ iij) 



> C vq 



(Al) 



where e = e + ji^sgne and the angular brackets imply 
averaging over the Fermi-line. To perform this averaging 
one can use the following identity: 



vq — te 



«2 



1 - 



r/(e,9). (A2) 



There are six possible configurations of the poles which 
give no n-zero contributions to the integral over ^ in 
Eq.(Al). Straightforward calculation yields the follow- 
ing expression for the current vertex 



7£(q;f^,w) 

where 

?7i ie,rt,uj) 



27rA^(0) q 



?7,(e,fi,w)/(e,,g) , (A3) 



1 + CJT SI 



— (die)e(uj~e)9(e-n) 
me V 



+ 0{-s)e{s - oj)e{n - e)y 



(A4) 



1 + UJT sgn £ V 



+ 9{-e)0ie ~ Lu)9{e - n)^ 



(A5) 



773(e, n,uj)= - r/4(£, n, uj) = (0{e)e{e - u)B{e - n) 

- 9{-e)e{uj - e)9{n ~ e)'^ (A6) 



and 



ei = 64 = 2e — rj + r ^ sgne, 
€2 — £3 — 2e — il — UJ + T^^ sgne. 



(A7) 
(A8) 



This presentation is convenient when w e cal culate the 
current vertex in the magnetic field (Sec. HE). 
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In the dirty limit Eg. ([A^) reduces to the local equa- 
tion, since in the limit r — > 



f{e,q)^Vq\ 
and thus the current vertex reads 

7(q) = cq = -47rA^(0)X'r^q. 



(A9) 



Here we keep only the terms which do not contradict to 
the thcta-fimctions in t he coopcrons, i.e. the third and 
fourth terms in Eq.(A3). Wc sec, that in the coordinate 
representation the vertex has the form: 

7(r) = -cV(5(r). 

This (5-functional behavior implies the locality of the cur- 
rent vertex in the dirty case. 

Let us also note, that if the external frequency cu is 
zero, the current vertex is easily connected with the 
particle-particle bubble for any r: 



1™ 7e(q;f^,w) 



9q ' 



(AlO) 



where 



n,(q; n)^N{Q)/^l d^g,{p)ga-e{ti- p] 

9{eie-n)) 



27r7V(0) 



'2£-17+isgne) 



(All) 
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